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Abstract 

We introduce the notion of parafermionic fields as the chiral fields which describe parti- 
cle excitations in two-dimensional conformal field theory, and argue that the parafermionic 
conformal dimensions can be determined using scale invariant scattering theory. Together 
with operator product arguments this may provide new information, in particular for non- 
rational conformal theories. We obtain in this way the field theoretical derivation of the 
critical exponents of the random cluster and 0(n) models, which in the limit of vanishing 
central charge yield percolation and self-avoiding walks. A simple derivation of the relation 
between S'-matrix and Lagrangian couplings of sine-Gordon model is also given. 



1 Introduction 



Field theory is the natural framework for the study of universal properties of statistical models 
exhibiting continuous phase transitions. The two-dimensional case remains exceptional to date 
because it allows for exact solutions. In particular, two-dimensional conformal field theory pQ 
yields a classification of renormalization group fixed points which can be associated to critical 
points of statistical models through symmetry arguments, with the determination of critical 
exponents as a main consequence. 

In two dimensions exact solvability is often possible already on the lattice [2] and many 
critical exponents provided by conformal field theory were first obtained within some discrete 
formulation of the given universality class. There are, however, examples for which the lattice 
derivation has been the only one available to date. This is the case of the two basic universality 
classes associated to S q (permutational) and 0(n) symmetry. The first one is represented on the 
lattice by the (/-state Potts model, which allows for a representation with q real corresponding 
to a model of random clusters [3], with ordinary percolation arising as q — > 1. The 0(n) model 
also allows for continuation to real values of n and gives self-avoiding walks in the limit n — > 

01- 

The critical exponents of the Potts and 0(n) models have been obtained within the frame- 
work of the lattice Coulomb gas (see [S] for a review). This consists of a number of mappings 
between lattice models leading to an identification with a system of interacting electric and mag- 
netic charges, which in turn can be related by renormalization group arguments to a Gaussian 
model for which exponents are computable. While this non-trivial procedure makes a field theo- 
retical derivation highly desirable, the continuation to real values of q or n leads to non-rational 
conformal field theories for which a canonical conformal description is not available. This is 
why, when the lattice results for the critical exponents were cast into the conformal language in 
[6j[7], the essential issue of the dependence on the parameters q and n remained as an external 
input inferred from the lattice results. 

In this paper we give a purely field theoretical derivation of the critical exponents of the S q 
and 0(n) models, independent of any lattice input. This is obtained as a consequence of a more 
general observation about the relation between the spin of the fields associated to the critical 
particle excitations (which in general is fractional, justifying the name of parafermionic fields) 
and scale invariant scattering amplitudes. Together with the fact that we are able to determine 
the parafermionic spin from operator product expansion and mutual locality arguments, and 
the scattering amplitudes from the enforcement of S q or 0(n) symmetry, the above mentioned 
relation leads to a relatively simple determination of the exponents. 

The paper is organized as follows. In the next section we recall the notion of mutual locality 
and introduce parafermions and their scattering properties. Section 3 is devoted to conformal 
field theory with central charge c < 1, in particular to characterize the critical lines with c < 1 
with the properties required for the purposes of this paper. The critical scattering theories with 
S q or 0(n) symmetry are then solved in section 4. Section 5 summarizes the results and contains 
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some additional comments. 



2 Fields and particles at criticality 

Consider a spin model of classical statistical mechanics which undergoes a continuous phase 
transition when a temperature parameter T takes a critical value T c . Within the field theory 
describing the universal properties of the scaling limit we call energy field e(x) the field con- 
jugated to the parameter r ~ T — T c ; we denote by x = (x\,X2) a point in two-dimensional 
Euclidean space. Typically the phase transition arises from the spontaneous breaking of an 
internal symmetry and the energy field is invariant under the action of the symmetry. The spin 
(or order parameter) field cr(x), on the other hand, transforms according to a representation of 
the symmetry and has in general several components (we drop indices for the time being). 

The fields e and a are examples of scaling fields. In two dimensions a scaling field $ is 
characterized by the conformal dimensions A$ and A$, which determine the scaling dimension 
Xj> = A$ + A<j> and the spin s<j> = A$ — A<j>. The field 3? picks up a phase e~ JS * e under a 
rotation by an angle 9 on the Euclidean plane. The scaling dimensions of e and a (which are 
scalar, i.e. spinless, fields) determine the correlation length and magnetization critical exponents 
as v = 1/2(1 — A e ) and /3 = 2A a v, respectively; the other canonical exponents in zero external 
field are obtained from the usual scaling relations. 

The operator product expansion (OPE) of two scaling fields can be written as 

= Y, C ^ zAk ~ Al ~ A ^ Ak ~ A% ~ Aj ^^(o) , (i) 

where Aj = A$ i; {$fe} is a basis of scaling fields, are structure constants, and z and z are 
related to the Euclidean coordinates as z = x\ + ixi and z = x\ — ixi. The fields <£j and &j 
are said to be mutally local if the correlation functions containing the product $j(a;)$j(0) do 
not possess a branch point at x = 0, i.e. are single valued under the continuation z — > e 2i7r z, 
z — > e~ 2l7r z. This is the case if the condition 

7ij = s k - Si- Sj eZ (2) 

holds for all <S?fc giving a non-zero contribution to the r.h.s. of ([I]). The fields having a direct 
physical meaning for the statistical model should be local with respet to the energy field e(x). 

It is often useful to look at one of the Euclidean dimensions as imaginary time. Rotational 
invariance of the two-dimensional Euclidean system then becomes relativistic invariance from 
the (l+l)-dimensional point of view. At criticality, the elementary excitations of a relativistic 
(l+l)-dimensional theory are massless particles with sign of the momentum that cannot be 
changed by a Lorentz transformation, i.e. they are either right movers or left movers. These 
particles are described by chiral (or holomorphic) fields ip, with A^, = 0, if right movers, and vp, 
with A,0 = 0, if left movers; one also has A^, = A^ = = —s^. In a scale invariant theory right 
and left degrees of freedom are decoupled, so that if we consider the scattering of a right mover 
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Figure 1 : Graphical illustration of equation ([3]) . 



with a left mover there is no dynamical phase shif10. Right-left scattering, however, involves 
exchanging the position on the line of two identical particles, so that scattering entails in general 
a statistical factor. Since in absence of dynamical interaction the passage from the initial to the 
final state can also be realized by it- rotations (see Fig. [[]), the right-left scattering (statistical) 
phase can be written as 

g — e -^( s i>- s i,) = e -2i7rA _ /g\ 

Notice that A^, integer (half integer) yields bosonic (fermionic) statistics; the spin, however, can 
take more general values and for this reason we refer to tp and $ as parafermionic fields. 

In general parafermionic excitations and the corresponding fields form multiplets transform- 
ing under some representation of the internal symmetry of the theory. We denote by *ff a the 
components of such multiplets, and by ^f a their anti-chiral counterparts. We can exploit the 
fact that in right-left scattering the two particles exchange their relative positions on the line to 
encode the most general process in the symbolic relation 

V a o * b = S C J fjo$ c , (4) 

where are scattering amplitudes and summation over repeated indices on the r.h.s. accounts 
for the fact that in general the parafermionic components may not be scaling fields. In such a 
case, however, it should be possible to find a linear combination which diagonalizes @ yielding 
a scattering phase and, using (J3j) , a definite conformal dimension. 

The scattering amplitudes are subject to the constraints of unitarity and crossing symmetry, 
which in the present case simply read 



5cd 
ab 



ad 



m, (5) 

; (6) 



the asterisk denotes complex conjugation and the bar over indices denotes antiparticles. 

We will illustrate in the next sections that A^, can also be determined from the requirement 
of mutual locality between tjj and e, and that comparison with ([3]) may provide non-trivial 
information. 



1 Scale invariance prevents a momentum dependence of the scattering amplitude. The only relativistic invariant 
is the center of mass energy and there is no way of building a dimensionless scaling variable. 
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3 Conformal field theory with c < 1 

We will be interested for our purposes in investigating some general properties of conformal field 
theories with central charge c < 1. A way of doing this is to consider the conformal field theory 
of a free scalar boson with action 

A= ± fd 2 x(d a <p) 2 (7) 



4vr _ 

and components of the stress tensor given by 

T(z) = -(d<p) 2 + iQd 2 ip, f(z) = -(8ip) 2 + iQd 2 ip, (8) 

with d and d denoting derivatives with respect to z and z, respectively. The terms proportional 
to Q introduce a modification of the stress tensor which does not change energy and momentum. 
They affect, however, the central charge of the theory and the conformal dimensions of a field 
which appear as coefficients in the OPE of the stress tensor with itself and with respectively 
PP. The central charge becomes [7] 

c = 1 - 6Q 2 . (9) 

Introducing the decomposition <p(x) = 4>{z) + 4>(z) allowed by the equation of motion dd<p = 0, 
the primary fields correpond to the exponentials 

V P:P = e 2i W>+P<« t (10) 

with conformal dimensions (A p ,Ap), and [7] 

A p = p(p-Q). (11) 

c=l. It will be useful to start with the case Q = 0, namely the usual Gaussian fixed point 
theory with central charge c = 1 and OPE 

where we omit for simplicity coordinate dependence and structure constants; the square brackets 
indicate that, besides the primary field, the r.h.s. contains its descendants, with conformal 
dimensions exceeding those of the primary by positive integers. 

A generic choice for the energy field in this theory is e = & + V—b _& oc cos 2b(p with 
A £ = b 2 . Using (fT2|) , (|2j) and (JTTj) . we obtain that the fields V Pt p are local with respect to e if 
A p ±b — Ap±b — A p + Ap = ±2b(p — p) is an integer, i.e. if 

TYl 

p-p=—, m£Z. (13) 

It becomes transparent in the fermionized (Thirring) version of the theory [51 [9] that the integer 
m is the U(l) fermionic charge. The fundamental parafermions will be the chiral fields local 
with respect to e with the lowest charge, i.e. = V±i/2&,o ; ^± = ^b,±i/26> with A^ ± = 1/46 2 . 
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On the other hand, the spin field components associated to the U(l) ~ 0(2) symmetry are 
selected picking up the scalar fields with m = ±1, i.e. a± = V ± i^ b ^ T i/ 4b with A CT± = 1/166 2 . 

We see in particular that at b 2 = 1/2 the parafermions become ordinary spin 1/2 fermions, 
so that this is the point where c = 1 results from two decoupled Ising models, each with c = 1/2, 
energy E{ and spin <j{ (i = 1, 2). We have A e = A £i = 1/2 and A a± = A CT1(J2 = 2A CTj = 1/8. 

c<l. When Q is real the central charge ([9]) is smaller than 1. It is useful to define the parameters 
/3 and t through the relations 

Q = {3- 1 -[3, /3 2 = t/(i + l), (14) 
so that, in particular, c = 1 — j^pyy- We also introduce the notations 

p*v = \[Q--n)p- 1 -o--m, as) 

*w(z) = V p ^,o , = Voj,^ ; (16) 

it follows from (jlip that the single non- vanishing conformal dimension of the fields (|16p is 

_ [(t + l^-ti/] 2 -! 
^~ 4i(i + l) • [U) 

It is a fundamental result of c < 1 conformal field theory that when the indices fi and v take 
positive integer values m and n the fields (fTS]) become degenerat^, and are characterized, in 
particular, by the degenerate-nondegenerate and degenerate-degenerate OPE's [T] 

m— 1 n— 1 

} , (18) 

fc=0 z=o 

min(mi,m2)-l mm(ni,n2)-l 

$mi,m ' ^m 2 ,n 2 = [^|mi-m 2 |+l+2fc>i-n 2 |+l+2j] ! (19) 

fc=0 Z=0 

similar relations hold for the fields ^p,^- 

The energy field e is consistently defined if the OPE e • e does not produce fields with 
scaling dimension smaller than X £ . It does not seem possible to satisfy this condition if e is non- 
degenerate, and it can be ckecked that the only degenerate fields which satisfy this condition, are 
relevant (in the renormalization group sense, namely have A < 1) and have positive dimension 
in the interval < c < 1 are $1,2, $2,1 an d $1,3. We will confine our discussion to theories 
which for c = 1/2 (t = 3) reduce to the Ising model with A £ = 1/2; then we focus on $2,1 and 
&13 which indeed satisfy this requirement. 

Let us start from e = $2,1^2,1 an d look for the corresponding parafermion ip{z). One can 
use (|18|) and ([2]) to check that there are no non-degenerate chiral fields local with respect to e. 



Rational values of (3 yield rational conformal theories that can be consistently built on a finite number of 
degenerate primaries and their descendants [T]. 
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On the other hand, for a degenerate chiral field (|19j) gives 

min(2,m) — 1 

S ■ $m,n = [*|m-2|+l+2Jfe,Al] . ( 2 0) 

k=0 

and ([2j) becomes 

A| m _ 2 |+i+2fc,n - A 2> i - A m , n = -_[2m-2A:(l + fc)-|m-2|(l + 2fc)-l] 

+-[l-n(l + 2fc-m+|m-2|)] G Z. (21) 

These conditions are satisfied for any /3 if m = 1 and n is odd; the requirement = 1/2 at 
t = 3 finally fixes n = 3, i.e. V = ^1,3- 

When looking for the spin field it is natural to require for generic c the condition 

e . a = a + ... (22) 

which holds in the Ising case. Since (fT8|) gives $2,1 ' ^fi,u = [^^-i,u] + [^+1,1/] i ([22]) is satisfied if 
a = ^1/2,0^1/2,0 an d we identify the fields and Q-^-y, which have the same dimension. The 
OPE (|22p also generates the scalar field ^3/2,0^3/2,0 which is naturally identified as a subleading 
spin field; more generally, iteration of the OPE with e produces the family of subleadingspin 

fields cjfc = $/fc+i/2,o^fc+i/2,o> k = 1, 2, On the other hand, e • e and iterations producqj the 

subleading thermal fields Ek = 3>fc+2 l^fc+2 l> fe = 1, 2, . . .. The structure of the OPE is invariant 
under the transformation — )• (— l) fc+1 £fc, including eq = e. 

We can use (fT8|) to see that the OPE ip ■ a = $1,3 • ^1/2,0^1/2,0 produces a field, let us call 
it fj,, with the same conformal dimensions of a. This in turn leads to a ■ fi = tp + • • • , meaning 
that a and pL cannot coincidj^l because they are mutually non-local. The existence of the field 
fi is not surprising once we recognize Ek — > (— l) fc+1 £fc as a duality transformation under which 
a — > fj,; such a duality is well known for the particular case of the Ising model. We see that 
for e = $2,1^2,1 the notion of parafermion we introduced matches the usual one |1CH lll| of field 
produced in the OPE of the spin field with its dual. 

We conclude this section with the case e = $1,3^1,3 . An analysis completely analogous to 
that of the previous case selects ip = $2,1- The OPE e ■ £ and iterations generate the subleading 
thermal fields e^ = 3>i 2fe+3^i 2fc+3 ; A; = 1,2,..., without that a dual symmetry emerges in this 
case. As for the spin field, this time the condition ()22[) is automatically satisfied and gives no 
constraint. 

3 The requirement that s is local with respect to itself implies that only scalar primaries and their descendants 
are actually produced in the OPE e ■ e. A similar prescription has to be adopted for the OPE's of the fields a^.. 
and £fc with s. 

4 This is possible beacause the conditions for degenerate fields only determine the dimensions of the fields 
appearing in the OPE. Fields with the same dimensions may differ if this is required by consistence with mutual 
locality. 
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4 Scattering and symmetries 



In the previous section we identified two critical lines of conformal field theory with central 
charge smaller than 1. The analysis was based on degenerate fields and their OPE's and does 
not seem to allow the identification of internal symmetries in any obvious way. On the other 
hand, internal symmetries play an essential role within the scattering framework described in 
section 2. We now apply the scattering formalism to the cases of 0(n) and S q symmetries and 
match the results with those of the previous section. 

4.1 0(n) symmetry and self-avoiding walks 

It is natural to start with a multiplet of neutral excitations labeled by an index a = 1, . . . , n. 
The relation takes the 0(n)-covariant form (Fig. [2]) 



*a o = 5 ab Si *c o * c + S 2 $ t o$ a + S 3 $„o$ 6 , (23) 

c=l 

The unitarity equations read 

S 2 S* 2 + S 3 S* 3 = 1 , (24) 
S 2 S% + S 3 S% = , (25) 
nS^l + S 2 Sf + SiS 2 * + S 3 Sf + SiSI = , (26) 

while crossing symmetry gives 

S l = S* 3 = pxe* , S 2 = S* 2 = p 2 , (27) 

with pi non-negative and p 2 real. The unitarity equations are then rewritten as 

pl + pl = l, (28) 
pip 2 cos(p = 0, (29) 
np\ + 2/31/32 cos ip + 2p1 cos 2(p = . (30) 

There are three ways of satisfying (f2T)j) . The first possibility is that cos</? = 0, so that (|3"0"|) fixes 
n = ±2; we will come back to this case in section 4.3. The second possibility is that pi = 0, but 
then we are left with S 2 = ±1 as the only non- vanishing amplitude, i.e. with n decoupled free 
bosons or fermions, an uninteresting case. Hence, the only non-trivial case with a continuous 
n-dependence is 

92 = , pi = 1 , n = — 2 cos 2tp . (31) 

It is natural to think of the particle trajectories as those of walks of n different colors on the 
Euclidean plane. The result S 2 = then amounts to a condition of self and mutual avoidance 
of the walks (Fig. [2]). This was first observed in [T2] in the study of the off-critical case. 
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Figure 2: Amplitudes Si, S2 and S3 of the 0(n)-invariant theory. Time runs vertically. 
It follows from fl23]) that 

n n 
^2^ a o^ a = S ^* a O* a , (32) 
a=l a=l 

S = nSi + S 2 + S 3 = -e 3 ^ , (33) 

where in the last equality we used the solution ([3T]) . Equation ([3]) then relates ip to A^, up to 
the 2nk ambiguity involved in the comparison of phases. 

The solution ([3Tj) covers the range n G (—2,2). The relations (f3Tj) and ([33]) consistently 
give n = 1 and the fermionic (Ising) value S = —1 when 93 = 27rfc/3, = 1,2 (mod 3). The 
0(n)-invariant critical line we are discussing should then be matched with one of the two critical 
theories with c < 1 that we identified in the previous section. Since no order-disorder duality is 
expected for 0(n) symmetry, the natural choice is the theory with e = $1,3^1,3 and ip = $2,1 ■ 
We know that A e = Ai^ = 1 at c = 1 and that for A e = 1 (i.e. for b 2 = 1) on the 0(2) line 
the dimension of the parafermion is Aw, = 1/4. With this further piece of information we fix 
93 = — ^-(A^, + 1/2) (mod27r); substituting A^ = A2,i and plugging into (f3TJ) we finally obtain 

n = 2 cos ~ , (34) 

where t was introduced in (|14p and parameterizes the central charge. 

We saw that the spin field has A CT = 1/166 2 along the 0(2) line with c = 1, and that the 
n = 2 endpoint of the solution (I3ip corresponds to b 2 = 1. Hence we conclude that for the latter 
theory A CT = 1/16 not only at c = 1/2 but also at c = 1; using these conditions to determine the 
indices (assumed ^-independent) in (|17p gives A CT = Ax/2,0 as ^ ne only positive solution within 
the range c G (0, 1). The iteration of the OPE with e then produces subleading spin fields 
with dimension A 1 / 2 ,±2fc; 

k = 1,2,.... 

4.2 S'g symmetry and cluster boundaries 

The natural particle basis for an S g -invariant theory below the critical temperature is that of the 
massive kinks corresponding to domain walls between the q different ordered phases produced by 
spontaneous symmetry breaking [13]. Such an interpretation of the excitations does not apply 
as it is to the critical case we want to investigate because at criticality the symmetry is unbroken 
and there is a single disordered phase. On the other hand, we can keep the same representation 
of the symmetry if we interpret the trajectories of the massless particles as the boundaries 
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Figure 3: Amplitudes So, Si, S2 and S3 of the S„-invariant theory. 



separating clusters of spins of different colors. We then denote by ty a p {a, j3 = 1, . . . , q; a ^ j3) 
the parafermion corresponding to the particle whose trajectory separates a cluster of color a 
from one of color (3. Permutational symmetry of the colors now gives to @ the form (Fig. [3j) 



m ai o ^ 7/3 = (1 - 5 a/3 ) 



S *<*<5 O V S p + Si $ Q7 o * 7/3 

S 2 ^ ^Sa + S3 *a 7 O 

5^7 



(35) 



where it is understood that the two indices of a parafermionic field have to be different. 
Crossing symmetry leads to the relations 



So — Sq = po 



Si = S 2 * = pe* , 



S3 — S3 — /?3 , 



(36) 



with po and P3 real and p non-negative. With these parameterizations the unitarity equations 
read 



(q - 3)p 2 + p 2 = l, 

(q - 4)po + 2poP cos ip = , 

(q - 2)p 2 + pl = l, 

(q - 3)p 2 + 2pp 3 cos ip = . 



(37) 
(38) 
(39) 
(40) 



It is straightforward to see that they produce a quadratic equation for p 2 with solutions p 2 ^ 
4 — q and p%\ = {q — ^)/{q 2 — 5g + 5). Choosing the first one leads to 



Po 



-1 



2 cos ip = — \/ 4: — q , 



P3 



(41) 



provided we use the condition ^3 = — 1 at q = 2 (Ising) to fix the overall sign of p^. It follows 
from the form of the solution (|^T|) that it holds for q 6 (0, 4). 

An interesting way to see that (I41j) is the solution we are interested in is to observe that a 
cluster boundary can branch into two boundaries at a point where three clusters meet, so that 
we expect the fusion relation 

• ^7/3 ~ > ( 42 ) 

with a structure constant C*^ which is color-independent by permutational symmetry. The 
requirement that the result of the scattering of $?ps with the r.h.s. equals that of the scattering 
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Figure 4: Diagrammatic equation which for different choices of external indices produces (|43H46|) . 



with the l.h.s. yields the equations (Fig. H]) 

51 = Si + {q- 3)S 2 S , So = S 3 S + Si + (q- 4)S 2 S , (43) 

S 3 = (q- 3)S 2 , S 2 = Sj + (q- 3)S S 1 , (44) 

5 2 = Si So +(q- 4)S 2 , Si = S S 2 + (q- 4)S 2 , (45) 
So = S 2 + S S 3 + (q - 4)S S l , So = Si So + S S 2 + (q- 5)S 2 , (46) 

whose unique non- vanishing solution is (|4ip . 

At q = 2 the only physical amplitude is S3 = —1, as required for the Ising model. The 

solution (|4ip then corresponds to a critical line spanning a neighborhood of c = 1/2, and it is 



naturally associated to the remaining theory with c < 1 of the previous section, namely the 
one with e = $2,1^2,1 and ip = $1,3; the order-disorder duality that we found in such a theory 
is known to hold for the two-dimensional Potts model. The maximal value of central charge, 
c = 1, should correspond to the maximal value of q, which we saw is 4. 
It follows from (f35|) that 

^«7 ° ^7" = S Yl ° * T Q ' ( 4? ) 

7 7 

S = S 3 + (q - 2)S 2 = e" 4 ^ , (48) 

with the last equality produced by the solution (|4ip . Comparison with ©, together with 
the conditions that = Ai ; 3 equals 1/2 and 1 for q = 2 and 4, respectively, determine 
ip = 7r(l + tjAi ; 3) (mod27r), which substituted in (jiTj) leads to the relation between q and the 
central charge 

v ^ = 2sin 7r / ( *~\ ) . (49) 
vy 2(t + l) V 7 

4.3 0(2) symmetry and sine-Gordon coupling 

We saw that (|28ti30|) admit a solution with n = 2. It can be parameterized as 

pi = sin a, p2 = cosa, ip = —7r/2, (50) 
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Symmetry 


c 






A CTfc 


0(2) 

0(n) n = 2cosf 
S q ^ = 2sm 2( ( m j 


1 

1 6 

1 t(t+i) 

1 6 

1 t(t+i) 


6 2 

Al,2fc+3 
Afc+2,1 


i 

W 
A 2 ,i 

Ai, 3 


l 

Al/2,±2fc 
Afc+1/2,0 



Table 1: Central charge and conformal dimensions along the three critical lines analyzed in the 
paper. A^ is the parafermionic dimension, while A £fe and A CTfc are the thermal and magnetic 
dimensions, respectively, with the leading dimensions (the only ones quoted for 0(2)) corre- 
sponding to k = and the subleading ones corresponding to k = 1,2,...; A^ )t , is given in 

(HZD. 



where we now allow for p\ negative. The phase in (j32|) becomes 

S = 25i + S 2 + S 3 = e~ ia , (51) 

and has to be equated to ^ with A^, which, as seen in the previous section, takes the value 
1/46 2 along the 0(2) critical line; this gives 

a = W- (52) 

The S-matrix of a massive 0(2)-invariant integrable theory was determined in p3] and found 
to depend on a non-negative parameter £. This S'-matrix corresponds to the sine-Gordon model 
obtained perturbing the 0(2)-invariant critical line with c = 1 by the field e = cos 2bip, and £ is 
a function of 6; the S*-matrix becomes minus the identity when £ = ir, which then corresponds 
to the free fermion point b 2 = 1/2. The elementary excitations are solitons and anti-solitons 
interpolated by the fields *$>± = ± i^2- As a consequence, the asymptotic high energy limit 
of the soliton-antisoliton amplitude, which turns out to be e 21 [H], should coincide with 
the massless amplitude (f5T|) . Using (f52|) and the condition at £ = it one obtains the relation 

(53) 

which was introduced in p3] through semiclassical and perturbative arguments. 



5 Conclusion 

The relations (|34p and (|49p determine the n- and g-dependence of the conformal dimensions 
identified in section 3 and at the end of section 4.1; the results are summarized in Table 1. The 
values of the leading and some subleading dimensions A £fc and A ak obtained within the lattice 
Coulomb gas framework [5] can be checked to coincide with those in Table 1. The latter also 
gives the information about the 0(2)-symmetric critical line. It was shown in section 4.3 how the 
formalism of this paper allows for a simple derivation of the relation (|53p between the coordinate 
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b along this line and the parameter £ entering the S- matrix of the massive 0(2)-invariant theory 
(sine-Gordon model). 

Other symmetries can be considered. For example, Z^ is implemented on a basis of parafermions 
ty a , with a = 1, 2, . . . , N — 1 denoting the Z^r charge and A^ Q = A^ N _ a by charge conjuga- 
tion. The diagonal solution S^f = 5^5 b S ab is subject to the unitarity and crossing constraints 
Sab = ^/S* ab = S* aN _ h . Requiring the fusion ^ a • ^ b ~ ^ a +b(modN) leads to 5 ac 5 bc = S a +b,c and 
to the result S ab = e 2inab / N , which in turn gives (i/27r) In S aa = -a 2 /N = A# (modi). The 
requirements that = 1/2 in the Ising case N = 2 and that for N = 3 A«j, 1 = A^, 2 takes the 
value 2/3 that we obtained above for S3 symmetry lead to A$ a = a(N — a)/N . This value is a 
starting point of the construction of Zjy-invariant conformal field theories in [11] . 

This work originates from a study concerning the space of fields of massive field theories [15] 
where, in particular, it was observed that parafermionic fields naturally appear in the analysis of 
high energy asymptotics. In the present paper, working directly at criticality, we have seen how 
particle excitations are described by chiral fields which in general have parafermionic nature, and 
how non-trivial results, in particular for non-rational conformal field theories, can be obtained 
studying parafermions in both the operator product and scattering settings. It is interesting to 
notice that in recent years discrete parafermionic observables were found to provide an intriguing 
link between criticality, holomorphicity and integrability on the lattice (see [16] for a review). 
The universal role of parafermions exhibited in this paper directly in the continuum may also 
be relevant for the interpretation of some of the lattice findings. 
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